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Introduction 

Grant NGR 33-016-174 was initiated at New York University 
(NYU) on 1 November 1971 for a nominal period of one year, and 
extended for a second year with an expected anniversary date of 
31 October 1973. Because of the closing of the NYU School of Engi- 

V. 

neering and Science on 31 August 1973, the grant has been prematurely 
terminated as of the latter date. By accelerating our efforts in anti- 
cipation of the NYU closing, we have been able to accomplish most of 
our research objectives. However, as noted below, two of our prin- 
cipal research results will not be ready for publication before the 
new termination date, although all the computations and most of the 
figures have been completed. This report does, nevertheless, pre- 
sent a brief account of all the research completed since the inception 
of the grant. 

The major purpose of the research project has been to carry 
out two main tasks in cooperation with the Goddard Institute for Space 
Studies (GISS). The first task was to provide assistance to GISS in 
its work on the development and application of numerical-dynamical 
models of the atmospheric global circulation. The second task has 
been to answer certain basic questions regarding the feasibility of 
long-range weather prediction by dynamical methods. The principal 
effort under the second task was devoted to studies of long term 
responses of the model atmosphere to anomalies in snow cover and 
sea surface temperature (SST). 
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Technical support to GISS, as part of the first task, was 
carried out primarily through conferences and visits to GISS by the 
principal investigator. Experimental outputs from several dynamical 
prediction models developed at GISS were analyzed in concert with 
the GISS staff to determine if the models faithfully simulated pro- 
perties of the real atmosphere. Most of this work was associated 
with the efforts at GISS to develop a global nine- level atmospheric 
model on the foundation of the original two-level model developed at 
UCLA by Mintz and Arakawa. No reports or publications by the NYU 
project staff resulted from this phase of the project, as our contri- / 
butions were mainly consultative and advisory. 

Most of the project work has been devoted to the second task, 
with primary emphasis on the response of the two-level Mintz - 
Arakawa general circulation model to pools of anamalously warm 
SST in the North and South Pacific Oceans. Three reports and two 
publications, plus a third probable publication, have resulted thus far 
from these efforts. Two additional reports, by research assistants, 

S. H. Chow and L. Lewis, are now being completed as a Ph. D. thesis 
and M. S. thesis respectively, and will be submitted for publication 
at a later date. In this final report, the results already reported and/ 
or published are presented only in the form of abstracts. However, 
a more extensive account of the work by Chow and Lewis is given 
here, as that work has not been previously reported. 

Most of the computations for this research project were carried 



out on the IBM 360/95 computer at GISS through the courtesy of 
Dr* .Robert Jastrow, director, and Dr. Milton Halem, and with the 
aid of Mr. John Liu. Some preliminary computations were per- 
formed on the UNIVAC 1108 computer at NYU. 



Reports and Publications 


The following technical reports have already been distributed. 

Jerome Spar, "Effects of Surface Anomalies on Certain Model 
Model Generated Meteorological Histories". NYU Geo- 
physical Sciences Laboratory Report No. GSL-TR-72-2, 
April 1972. 54 pages. 

Jerome Spar, "Transequatorial Effects of Sea Surface Tem- 
perature Anomalies in a Global General Circulation Model’ 
NYU Geophysical Sciences Laboratory Report No. GSL- 
TR-72-7, September 1972. 28 pages. 

Jerome Spar, "Supplementary Notes on Sea Surface Tempera- 
ture Anomalies and Model-Generated Meteorological 
Histories". NYU Geophysical Sciences Laboratory Re- 
port No. GSL-TR-72-9, December 1972. 24 pages. 

The following paper has been published. 

Jerome Spar, ’Some Effects of Surface Anomalies in a Global 
General Circulation Model", Monthly Weather Review , 

Vol. 101, No. 2, February 1973, pp, 91-100. 

The following paper has been accepted for publication in the July (or 
August) 1973 issue of the Monthly Weather Review . 

Jerome Spar, "Transequatorial Effects of Sea Surface Tem- 
perature Anomalies in a Global General Circulation 


Model". 



The following paper has been submitted for publication to the Mo nthly 

j 

Weather Review , and is now in the editorial review process. 

Jerome Spar, "Supplementary Notes on Sea Surface Tempera- 
ture Anomalies and Model-Generated Meteorological 
Histories". 

A report on the project research was presented by Jerome Spar 
in the form of an invited paper at a conference on "Climatic Changes 
on Time Scales from a Month to Millenia" at the.Scripps Institution of 
Oceanography in La Jolla, Calif, on 15-17 November 1972. This paper, 
titled "Response of a Model Atmosphere to Monthly and Seasonal Sea 
Surface Temperature Anomalies''^ is described briefly in the review 
of the conference published in the Bulletin of the American Meteorological 
Society, Vol. 54, No<> 5 , May 1973, pp. 425-432. 



Comput ed Resp onses to Surface Anomalies in a Global Atmospheric Model 
(Abstracts of previously issued reports and publications) 

The Mintz-Arakawa 2-level general circulation model was 
used in a series of experiments to compute the response of the model 
atmosphere to (1) a positive sea surface temperature (SST) anomaly 
in the North Pacific Ocean in summer and in winter, (2) an identical 
anomaly in the South Pacific Ocean in the southern hemisphere winter 
season, and (3) anomalous northward and southward displacements 
of the northern hemisphere snow line over the continents. In each case 
computations were carried out for 90 "forecast" days. Results were 
studied in terms of the differences between anomaly and control histories. 
Time series of certain regional response indices, including area-aver- 
age sea level pressure and 600 mb circulation indices, as well as 30- 
day mean sea level pressure maps, were used in the analysis. 

The experiments showed significant interhemispheric effects 
after about one month, phase shifts of 1-2 weeks in major cyclonic 
developments, stronger reactions to sea temperature anomalies in 
winter than in summer, and marked influence of the snow line on the 
winter monsoonal pressure difference between the continents and the 
North Atlantic Ocean. 

The global response of the atmosphere, as simulated by the 
Mintz-A rakawa model, to a persistent, anomalous pool of warm sea 
surface temperatures in the extratropical Pacific Ocean was also 
examined in terms of the meridional pole-to-pole profile of the zonally- 


averaged 600 mb surface for periods up to 90 days. Following an 
initial hydrostatic inflation of the isobaric surface in the latitude of 
the warm pool, effects spread poleward within the hemisphere, then 

began to appear after about two to three weeks in high latitudes of the 

\ 

opposite hemisphere, but with little or no response in the tropics. 
The same sea temperature anomaly field generated a stronger re- 
sponse in winter than in summer, and a very different reaction when 
located in the Southern Hemisphere than when in the Northern Hemi- 
sphere. After a month of thermal forcing, the response to the SST 
anomaly was at least as large in the opposite hemisphere as in the 
hemisphere of the anomaly. The winter hemisphere responded more 
rapidly to the SST anomaly in the opposite hemisphere than did the 
summer hemisphere.. Vacillation between low and high meridional 

wave number patterns was observed in the computed reactions to the 
warm pool. 

In seasonal computations, the Mintz-Arakawa model was 
found to be sensitive to a minor alteration in the computational pro- 
gram. Effects of the program change on monthly mean sea level 
pressure fields were small in the first month, but large in the second 
and third months, although the meteorological histories generated by 
both the original and modified programs were equally credible. 

The inherited effects of a "transient’ 1 SST anomaly (lasting 
only one month) on the computed monthly mean sea level pressure 
fields over the period of a season were about as large in absolute 
magnitude as those generated in the model by a "persistent" SST 
anomaly lasting three months. 


The effects of the transient SST anomaly in the North Pacific 
Ocean on monthly mean temperature and precipitation in the eastern 
United States were large enough to produce changes corresponding to 
one or two of the class intervals used by the National Weather Service 
for these '’predicted ' 1 weather elements. The model-generated pre- 
cipitation in the equatorial region was also found to be sensitive to the 
SST anomaly in the North Pacific. 
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Efjects of a Transient Sea Surface Temperature Anomaly on the 
Energetics of the Mintz-Arakawa Model Atmosphere 

v 

(Shu-Hsien Chow) 

Some experiments concerning the effects of sea surface tem- 
perature (SST) anomalies have been carried out at the Goddard In- 
stitute. for Space Studies with the two-level Mintz-Arakawa global 
general circulation model (Gates, et al. , 1971)* Certain hypothetical 
anomaly patterns were superimposed on the sea surface temperature 
field of the Pacific Ocean in those experiments with the anomaly pattern 
either persistent over a season or over one month* The response of 
the model atmosphere to the SST anomalies were discussed in terms 
of the monthly mean sea level pressure patterns and the regional in- 
dices. In order to gain more insight into the effects of the anomalies, 
we have also investigated them from the energetics point of view. 

The energetics of the atmosphere have been investigated in- 
tensively during the past few years. It is generally confirmed that the 
major energy source for the atmosphere is the generation of zonal 
available potential energy by diabatic processes (e. g. , Oort, 1964, 
Dutton and Johnson, 1967). Most of the zonal available potential ' 
energy, which is generated by low -latitude heating and high- latitude 
cooling, is converted into eddy available potential energy ( Wiin-Nielsen, 
Brown and Drake, 1963, 1964), while only a small energy exchange 
takes place between the zonal available potential energy and the zonal 
kinetic energy (Wiin-Nielsen, 1959, Saltzman and Fleisher, I960, 

1961). 



All the observational studies indicate that energy is trans- 
formed from eddy available potential energy to eddy kinetic energy 
through the rising of warm air and the sinking of cold air around a 
latitude circle (Wiin-Nielsen, 1959, Saltzman .and Fleisher, I960, 
1961, Kruegejj Winston and Haines, 1965). A large fraction of the 
eddy kinetic energy is dissipated by frictional processes, while the 
remaining small part is converted into zonal kinetic energy for main- 
taining the zonal wind system (Wiin-Nielsen, Brown and Drake, 1963, 
1964). 


The basic nature of the energetics of the atmosphere may be 
changed in unusual circumstances. For example, diagnostic studies 
the kinetic energy exchange between zonal flow and large-scale 
eddies for the troposphere north of 20°N have indicated the existence 
of an unusual type of circulation in January 1 9 63 (Wiin-Nielsen, Brown 
and Drake, 1964, Murakami and Tomatsu, 1965, Brown, 1967). Diag- 
nostic calculations for that month indicate that kinetic energy was 
transformed from zonal flow to large-scale eddies during this winter 


period, whereas this exchange mechanism is 
direction as mentioned before* 


usually in the opposite 


The abnormality of the circulation in January 1963 has been 
Studied particularly by Namias (1963) and O'Connor (1963). Accordi 


to Namias (1963), positive 
north Pacific Ocean were 
anc! early 1963. The case 


SST anomalies in the central and eastern 
quite persistent during the last half of 1962 
of January 1963 seems to suggest that SST 



anomalies may produce significant effects on the energetics of the 
atmosphere. Following the hints provided by these synoptic and diag- 
nostic studies, we have undertaken to calculate the effects of an SST 
anomaly on the energetics of the model atmosphere over a period of a 
season. For these model diagnostic calculations we have used the 
case of the transient (one month) SST anomaly in the North Pacific 
Ocean* 

The equations for studies of the energetics of the whole global 
atmosphere in the wave-number domain were first treated systemati- 
cally by Saltzman (1957) using one -dimensional Fourier analyses 
around latitude circles. The method of zonal harmonic analysis is also 
adopted in this study. The energy equations for the 'nth component 
of eddies and for the zonal flow, per unit mass over a latitude belt, 
may be written as 

(K(n)) t = F(K(n)) + < K(n) > + < K(m), K(n) > 

+ < P(n), K(n) > + W(n) - D(n), (1) 

(P(n)) = F(P(n)) + < P , P(n) > + < P(m), P(n) > 

j L Z 

- < P(n), K(n) > + G(n), (2) 

(K A = F < K Z > ' < K z* K E > +< P z - K z > + W z - D z , (3) 

< P z>t =F ( P z >- <r V P E > - <P z> K z >+G z 


(4) 



where ( ) t means the partial derivative with respect to time, K(n) 
is the eddy kinetic energy of the nth component, P(n) is the eddy 
available potential energy of the nth component, is the zonal ki- 
netic energy, is the zonal available potential energy,. F(K(n)) is 
the flux convergence of eddy kinetic energy of the nth component, 
F(P(n)) is the flux convergence of eddy available potential energy of 
the nth component, F(K z ) is the flux convergence of zonal kinetic 
energy, F(P z ) is the flux convergence of zonal available potential 
z > is the kinetic energy transfer from zonal flow 

to eddy of the nth component, < P z , P{n) > is the available potential 

energy transfer from zonal flow to eddy of the nth component, < K , K 

z E 

is the kinetic energy transfer from zonal flow to all the eddies, 

< P z* P E > is the available potential energy transfer from zonal flow 
to all the eddies, < K(m), K(n) > is the kinetic energy-transfer from 
eddies of other components to the nth component through nonlinear 
wave interactions, < P(m), P(n) > is the available potential energy 
transfer from eddies of other components to the nth component through 
nonlinear wave interactions, < P(n), K(n) > is the conversion from 
available potential energy to kinetic energy for the eddy of the nth com- 
ponent, < P 2 , K z > is the conversion from available potential energy 
to kinetic energy for the zonal flow, W(n) is the energy exchange with 
the surroundings through the eddy pressure interaction of the nth 
component, W z is the energy exchange with the surroundings through 
the zonal pressure interaction, D(n) is the frictional dissipation of 
the eddy kinetic energy of the nth component, is the frictional dis- 
sipation of zonal kinetic energy, G(n) is the generation of eddy avail- 
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able potential energy of the nth component by diabatic heating and 
is the generation of zonal available potential energy by diabatic heating. 

Note that by integrating Eqs. (1) - (4) over the whole global 
atmosphere, the resulting equations can be used to estimate the atmo- 
spheric energy cycle in the wave- number domain. Further, by sum- 
ming up over all the waves, the atmospheric energy cycle in the space 
domain can also be estimated. 

As mentioned above, the transient (one month) SST anomaly 

V 

experiment carried out by Spar (1972) was the subject of this study. 
Because the two-level Mintz-Arakawa general circulation model is 
designed to simulate the troposphere, only the energetics of the tropo- 
sphere are discussed. In the calculation of energetics, Lorenz's 
(1955) definition of available potential energy is adopted. Therefore 
p-coordinates are used instead of a -coordinates, (a = (p - p,j,)/(pg - Prj,), 
where p^ is the pressure of the tropopause and is taken to be 200 mb, 
and p Q is the pressure of the earth's surface). The troposphere is 
assumed to be bounded between 200 mb and 1000 mb, and equally 
divided into two layers with the 600 mb surface as the interface be- 
tween them. All of the variables at the p- surface are interpolated 
from those at the cr- surfaces. 

Velocities, frictional forces and diabatic heatings are inter- 
polated linearly in a- space. For the interpolations of temperatures 
and temperature lapse rates, it is assumed that potential tempera- 
ture is linear in p space, which is also adopted in the two- level 



Mintz-Arakawa general circulation model (p is pressure and k is 
the Poisson constant, 0. 288.) For the vertical p-velocities, a para- 
bolic profile in a -space is assumed* 

The energy transformation terms are calculated from the 
Fourier representations of variables defined at isobaric surfaces, and 
the harmonics are computed up to wave number 15, based on informa- 
tion at 72 points around a latitude circle. 

Two aspects of the problem under investigation are (1) how 
well does the model atmosphere simulate the atmospheric transforma- 
tion processes and energy contents, and (2) what are the effects of a 
transient SST anomaly on the energetics of the atmosphere. For the 
first problem, the energetics of the control run, which was designed 
to simulate the real atmosphere, are compared with observations. 

For the second problem, the energetics of the anomaly run are com- 
pared with the control run. 

Table 1 shows the model energy budget in the spatial domain 
for the winter control run for the region north of 20° N averaged over 
90 days, together with observational results for the same region from 
Krueger et al. (1965) and Saltzman (1970). The energy values of 
Krueger et al. (1965) are for an average of 5 winter (December - 
February) seasons, and are based on data from 850-500 mb, extra- 
polated to include the entire troposphere. Saltzman 's estimates are 
for the s Lx- mo nth "cold season", October - March, and are based on 
a variety of sources. 

The zonal available potential energy, P , of the model con- 
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Table 1 


Energy budget lor 90-day winter control run averaged 
over area north of latitude 20° No Observational values 
from Krueger et al. (1965) and Saltzman (1970) are shown 
for comparison. 

Model Krueger et Saltzman 


A, Energy contents 


p . 

670 

z 


K . 

130 

z 

P E . 

104 

K „ 

E. 

79 


(1 0 4 joules m ^ ) 


3 50 
152 
119 
134 


355 

100 

113 

143 


Energy Conve 

r sions 

< P , 

z 

K 2 

> . 

- 85 

<P Z' 

P E 

> o 

348 

<P E« 

k e 

> . 

589 

<K E’ 

K 

z 

> . 

51 


(10 ^ watt m ^) 


- 83 
513 
371 


0 

525 

303 

48 


Co 


Generations, Dissipations, and Flux Convergences (1 0" ^ watt m'^) 


G . 

z 

387 

430 

■ unvei wait m J 

300 

g e- 

| 62 

- 142 

- 159 

W o 

z 

247 

i 



E 

i 

! 1 



D . 

z 

268 

i 

j 

48 

D E- 

3 74 1 

i 


255 

F(P z ). 

154 ! 



F( p £ ). 

- i 



F( Kz ). 

4 



f(k e ). j 

11 
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troj J'un appears to be too large compared with the "observed" values. 
One reason for this discrepancy is that, in one sense, the domain 
considered for the model computations is different from that of the 
observational diagnostics* Squares of the temperature deviations 
from the global mean temperature were used to compute in the 
model computations, while deviations from the limited area mean 
temperature (north of 20° N) were used for the observational energetics 
Furthermore, the model produces excessively high temperatures in 
the upper troposphere, in the tropics, which also exaggerates P . 

z 

The zonal kinetic energy, , of the model is in good agree- 
ment with observations, and so is the eddy potential energy, P . 

E 

On the other hand, the eddy kinetic energy, K^, computed from the 

model is too low, probably as a result of the low horizontal resolu- 
tion of the model. 

The term < P^, > indicates a conversion from zonal kinetic 

to zonal potential energy in the model, in good agreement with the data 
of Krueger et al. (1965). The negative value of this conversion re- 
flects the dominance of the indirect Ferrel cell north of latitude 20” N. 
On the other hand, the conversion, < P £j K £ > , from eddy potential 
to eddy kinetic, appears to be too large in the model as compared with 
'bbservations". However, this conversion is calculated from a com- 
puted vertical velocity, which in turn is determined dry adiabatically 
in the "observational” diagnostics. Thus, the model, with its moist 
convection, may actually give more realistic values than the "observa- 


tions". 



The conversion < Kg, K ? > , from eddy to zonal kinetic 
energy, shows good agreement between the control case model and 
observations, but the term < P^ , Pg > , representing conversion 

from zonal to eddy potential energy, seems somewhat too low in the 
model. 

The generation of zonal available potential energy, G.. , by 

z 

diabatic heating in the model agrees well with the values estimated 
in the diagnostic studies, but the same cannot be said of the genera- 
tion of eddy potential energy, G^_. , for which there is a disagree- 
ment in sign. In view of the fact that G and G_ are. computed as 
residuals, rather than directly, in the observational diagnostic 
studies, they cannot be viewed as reliable standards against which 
to judge the model. 

The total kinetic energy dissipation, (D + D ), in the model 

Z 

- 2 — 

is 6.42 watts m compared with Saltzman^s residual estimate of 3. 03 

-2 

watts m . However, Kung (1967) has estimated a value of 4. 94 
- 2 

watts m for the winter dissipation north of 20° N, which is in better 
agieement with the model. It does appear likely, however, that D 

z 

is somewhat too high in the model. 

In general, the model control run exhibits a reasonable and 
realistic energy budget in the spatial domain that is well within the 
error limits of the diagnostic studies. 

The effects of the transient North Pacific SST anomaly on 
the gross energetics of the model in the spatial domain can be ex- 
amined in terms of the conversions and generations listed in Table 2, 
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T able 2* Model computations of energy generations and conversions 
for each winter month of control and anomaly runs in three 
zonal sections: 28°N - 60°N ("NH"), 20°N - 20°S ("Tropics"), 

and 28® S - 60°S ("SH"). 



Transient North Pacific 

SST anomaly. 

Units: 10 

watts m 

2 

« 



"NH" 

"TROPICS" 

i 

’SH" 


Control 

Anomaly 

Control 

A nomaly 

Control 

A nomaly 

Month No* 

1 

1 

1 

1 

1 

1 

G 

z 

177 

151 

932 

945 

118 

126 


78 

105 

840 

86 5 

52 

46 

<P , K > 
z z 

- 335 

- 307 

1088 

1095 

- 101 

- 91 

<P z* P E > 

555 

571 

- 162 

- 188 

99 

108 


1014 

1058 

813 

863 

66 5 

644 


79 

91 

21 

13 

47 

46 


Month No* 

2 

2 

2 

' 2 

2 

2 

G 

z 

291 ' 

323 

727 

713 

210 

197 

g e 

88 

0 

1407 

1327 

31 

13 

<P ,K > 
z z 

-100 

64 

1087 

1134 

- 95 

-138 

<p z > P E > 

512 

425 

- 356 

- 410 

78 

- 38 

<p P k e > 

852 

600 

1203 

1207 

556 

573 


11 

-8 

-14 

-30 

36 

59 
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T able 2 . (Cont'd, ) 



t 

"NH" 

'TROPICS” 


"SH 

r i 


Control 

A nomaly 

Control 

A nomaly 

Control 

A nomaly 

Month No 

. 3 

3 

3 

3 

3 

3 

G 

z 

231 

176 

861 

533 

185 

172 

g e 

84 

38 

1025 

167 

34 

28 

<P ,K > 
z z 

-112 

-191 

1080 

815 

-149 

-163 

<P z- P E > 

460 

363 

- 248 

- 30 

129 

87 


690 

729 

906 

155 

628 

561 


19 

65 

- 15 

10 

51 

82 


For each month of the 


experiment, the two generation values, G and 

° z 

Gg » anc * the four conversion values are given for both the anomaly and 
control run for three regions of the globe: 28°N - 60°N ("NH"), 20°N-20°S 
( Tropics ), and 28 P S - 60 ®S ("SH"). The units are the same as in Table 1. 

During the first month, the effects of the SST anomaly on the. re- 
gional generations and conversions of energy are quite small, partic- 
ularly in the tropics and Southern Hemisphere. In the Northern Hemi- 
sphere, the SST anomaly in the North Pacific appears to have caused 
a slight reduction in the diabatic generation of zonal potential energy and 
an increase in the generation of eddy potential energy, as might have been 

expected, together with slightly increased conversions from P to P 

z E 

P E to K E * and K E to K z - ^ The last result is contrarary to- the anomalous 
kinetic energy conversion of January 1963 }. The slight decrease in the 
magnitude of <P^,K^> reflects a weaker indirect Ferrel cell as one 
result of the SST warming. 

The effect of the transient North Pacific SST anomaly on the ' 
energetics is more apparent in the second month, and again mainly in 
the Northern Hemisphere. It must be recalled that the SST anomaly 
in this experiment was inserted as a step function, beginning as a dis- 
continuous SST warming on the first day of the first month and ending 
as a discontinuous SST cooling of the same magnitude on the first day 
of the second month. Thus, in the second month there are delayed 
carry over effects of the first month's warm pool, as well as effects 
of the SST cooling that returned the sea temperature field to the cli- 



;vuT,to.logical mean annual value u s u cl in the control run. Hence, a sim- 
ple interpretation of the resulting energetics is not possible. Never- 
theless, there are some noteworthy energy effects, particularly in the 
Northern Hemisphere winch are indicative of the disturbing effect of 
the transient warm pool. For example, in the anomaly run in the Northern 
Hemisphere the eddy generation term, Ci , has fallen to zero, the mo- 

It/ 

ridional generation term, < P , K > has been reversed in sign, indicating 
a positive, or direct, energy conversion cell in place of the negative, 
cr indirect, Ferrel cell, and the <K rr ,K > conversion has also changed 
sign, indicating an anomalous conversion from zonal to eddy kinetic energy. 
Outside the Northern Hemisphere, the most noteworthy effects of the 
anomaly in the second month are the reversal in sign of <P , P > in the 
Southern Hemisphere and the intensification of the indirect energy con- 
version process in the Southern Hemisphere represented by the increased 
negative value of <P^, K^>. These results for the second month seem 
to indicate that transient S5T anomalies can indeed have a disturbing cf 
enect on the global energetics of the atmosphere after a period of about 
cr.e month* 

The third month, as indicated in Table 2, is characterized by 
a return of the global energetics to the "normal 11 state represented by 
the control run (but with an important exception in the "Tropics"). 

In the Northern and Southern Hemispheres poleward of latitudes 28° 
(approximately the latitudes of the subtropical high pressure belts), 
the generation and conversion values for the anomaly run are not notably 



different from those for the control run, now that the soa temperature 
fields are the same. This does not mean that the extratropical synoptic 
meteorological patterns and structures are the same for the two runs. 

But it does indicate that the gross regional energetics in the spatial 
domain are once again similar. 

On the other hand, the energetics in the '’Tropics” are seriously 
distorted in the third month in the anomaly run compared with the con- 
trol run. Especially noteworthy are the order of magnitude decreases 
of G^, < > » and < P^» K > * There is good reason to believe that 
these latter results reflect a defect in the model, namely in the para- 
meterization of convection. In the tropics, the precipitation and latent •: 
heat release are due primarily to convection, which is very sensitive 
to the model's thermal structure. For some reason which has not yet 
been completely diagnos ed, the precipitation in the third month in the 
equatorial region was severely altered in the anomaly run compared with 
the control run. The convective precipitation has a dominating influence 
on the tropical energetics, with the results shown in Table 2. As of 
this writing we are of the opinion that this tropical effect is a com- 
putational artifact resulting from a. model defect, and not a valid reflec- 
tion of nature. 

A preliminary analysis of the model energetics in the wave 
number domain leads to the certain tentative conclusions regarding both 
the model's simulation of the spectral energetics of the real atmosphere 
and the effects of the SST anomaly on the spectral energetics. 

Except for the eddy generation terms, G(n), the model energy 
flow diagram in the wave number domain shows good qualitative agree- 
ment with the observational results of Saltzman (l‘>70). The available 



potential energy generated by diabatic heating takes place mainly in the 
zonal form, and flows to all the eddies mainly through eddy sensible 
heat transfer, with the dominant transfer in the larger eddies (smaller 
wave numbers). Available potential energy is further transferred from 
longer waves to intermediate and short waves through the non-linear 
wave interactions. All waves convert their potential energy to kinetic 
energy through baroclinic processes, represented by <P(n}, K(n)>, and 
supply kinetic energy to the zonal flow through the barotropic process 
represented by <K z> K(n)>. Intermediate waves lose kinetic energy, 
and most longer waves and short waves gain kinetic energy through the 
wave interactions, <K(m), K(n)>. However, the eddy energies, P(n) 
and K(n) are generally smaller in the model (except for n = 1 and 2) 
than in nature, as are the conversions from zonal to eddy form and the 
transfers of kinetic energy by wave interactions. Some of these short- 
comings of the model are undoubtedly due to the coarse horizontal re- 
solution, as shown by Manabe, et al. (1970). 

Maximum energy conversion, <P(n), K(n)>, in the model is 
found in the very long waves (near n = 2), whereas some theoretical 
and "observational" studies indicate a maximum in the cyclone scale, 
around n =6. However, those studies are based on dry adiabatic models 
or on diagnostic calculations for which the Vertical velocities were com- 
puted dry adiabatically. In the work of Haltiner (1967), Murakami 
and Tomatsu (1965), and Manabe, et al. (1970) it is shown that, when 
sensible heat transfer and other diabatic effects are included, maximum 
conversion does take place in the very long waves, as in the model. 

The fact that the model conversions at all wave numbers are larger than 
those reported by Saltzman is probably also due to the use of adiabatic 



vertical motions in the "observational" diagnostics. 

The effects of the SST anomaly on the spectral energetics are 
even more complex than on the gross energetics in the spatial domain. 

In this brief summary, we will confine our attention to effects in the 
Northern Hemisphere. 

In the first month, the SST anomaly is found to increase the 
eddy generation in the northern middle latitudes, mainly in waves 1 and 
2, due to the increase of eddy generation by heating and condensation 
and the reduction of eddy destruction by sensible heat exchange. The 
energy conversions, <P(n), K(n)> , in the intermediate- scale cyclone 
waves are found to increase in the anomaly run, probably because of 
increased baroclinic instability resulting fom the increased sea surface 
temperature gradient. Both eddy kinetic energy and eddy available 
potential energy are also found to increase at intermediate wave numbers 
in the anomaly run, and the kinetic energy transferred from eddy to zonal 
is found to increase as well. The kinetic energies transferred from eddy 
to zonal at waves 2 and 3 are found to decrease in the anomaly run. 

In the second month, removal of the SST anomaly from the 
anomaly run produces cooling effects. As expected, this decreases the 
heating and condensation, and reduces the eddy generation. The results 
show that the reduction occurs at all wave numbers with the major re- 
duction in the longer waves. The energy conversion is also found to 
decrease over almost all waves, with the major decrease in the longer 



waves. This is probably due to the decrease in vertical motion re- 
sulting from the decreased heating and condensation. Eddy kinetic and 
available potential energy are found to decrease at almost all the wave 
numbers in the anomaly run. 

In the third month, the combined effect of heating of the first 
month and cooling of the second month in the northern middle latitudes 
is found to result in increased eddy activity in the long and intermediate 
waves from n = 2 to n = 7. This can be seen in the increase of eddy 
kinetic energy, eddy available potential energy, eddy energy conver- 
sion, and the kinetic energy transferred from eddy to zonal at these 
wave numbers. In general, the spe ctral energetics in the Northern 
Hemisphere are as much disturbed in the third month as in the second 
month by the SST anomaly, unlike the gross energetics in the spatial 


domain. 
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The Angular Momentum Budget of the Mi ntz-Arakawa 

Modeli and Effects of a Transient SST Anomaly 

(L..J. Lewis) 

Asa further diagnostic test of the two-level model, and to 
gain additonal insight into the dynamical effects of a transient North 
Pacific SST anomaly, the zonal absolute angular momentum budget of the 
model atmosphere was computed for the 90 -day winter simulations with 

and without the anomalous warm pool. 

Due to surface friction (as well as pressure differences across 
mountains), there is a transfer of westerly (positive) angular momentum 

-®S 

from the earth to the atmosphere in the belts of surface easterly winds 
(mainly the trade wind region), and from the atmosphere to the earth in 
the mid -latitude belts of surface westerlies. To balance these additions 
and subtractions, there must be meridional transports of zonal angular 
momentum, particularly from the trades to the westerlies. In this study 
we have computed, as a function of latitude, each of the various contri- 
butions to the local time- variation of zonal angular momentum in toroidal 
atmospheric rings, integrated over the total depth of the model tropo- 
sphere. 

The terms in the angular momentum equation include (a) con- 
vergences of meridional fluxes of relative angular momentum by mean 
meridional motions and by eddies, (b) convergences of meridional fluxes 
of earth momentum by mean motions and by eddies, (c) convergences 
of vertical fluxes of relative and earth momentum by mean motions and 



by eddies, (d) momentum generation or dissipation by zonal pressure 
differences across mountains, and (e) momentum generation or dissipation 
by surface friction. 

Preliminary results of the calculations, based on the analyses 
completed thus far, indicate the following tentative conclusions. 

(1) The 90-day average meridional profile of the meridional 
momentum flux by mean motions in the Northern Hemisphere in the model 
control run is in good qualitative agreement with observational results. 
However, the maximum poleward flux in the model, at 16° - 20° N, is 
about 50% higher than in the mean winter atmosphere, apparently due 

to an excessively strong Hadley circulation in the model control run. 

(2) Meridional fluxes by eddies in the model are of the same 
magnitude, generally, as in nature. However, the profiles indicate that 
the width of the westerlies is narrower in the model control run than in 
the mean winter atmosphere in the Northern Hemisphere. There are 
marked changes in the eddy fluxes from month to month in the model 
control run, indicating changes in regime from low to high zonal index. 

(3) The mountain term in the model control run does not agree 
with that computed from the observed mean winter Northern Hemisphere 
This discrepancy appears to be due to the presence of a very strong 
Himalayan anticyclone in the model. 

(4) There is generally good qualitative agreement between 
the model and "obs erved" profile s of the surface stresses, and, hence, 

of the frictional torque terms. However, in the model, the strong Hadley 



cell if associated with stronger than average trade winds, and, there- 
fore, higher values for the surface stresses and frictional torques. 

All the comparisons with 'bbservation" mentioned above are 
for the Northern Hemisphere, as there is a paucity of data in the Southern 
Hemisphere. However, in the comparison of the model control and anom- 
aly runs, we are able to examine world-wide effects. Although the anal- 
ysis of the anomaly run is still not complete, the following effects of the 
SST anomaly are indicated. 

(5) The meridional fluxes by eddies are slightly increased and 
the fluxes by mean motions slightly decreased in the anomaly case re- 
lative to the control run, when averaged over the whole 90-day period. 

This effect is greatest, however, in the third month when, in the North- 
ern Hemisphere, the mean flux is reduced about 25% and the eddy flux about 
50% below the control run, bringing the anomaly profile into closer agree- 
ment with ’’observations". The surface stress profile also shows better 
agreement with observations in the anomaly run than in the control run.- 
These results appear to be due to a reduction in the strength of the 
Hadley cell in the case of the transient North Pacific warm pool com- 
pared with the control. The reasons for this behavior are not yet under- 
stood. 

A complete report on the angular momentum study will be pre- 
sented in the M. S a thesis by L. J. Lewis, and will be published or dis- 
tributed at a later time. 



